Are the low-momentum gluon correlations semiclassically determined? 
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We argue that low-energy gluodynamics can be explained in terms of semi-classical Yang-Mills 
solutions by demonstrating that lattice gluon correlation functions fit to instanton liquid predictions 
for low energies and, after cooling, in the whole range. 
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I. INTRODUCTION II. INSTANTON LIQUID PICTURE 



Semi-classical methods exploiting the non-dispersive 
solutions of non- linear classical field equations, named 
solitons, are widely used in physics. In particular, the 
time-dependent four-dimension solutions of the Yang- 
Mills equations (instantons, mcrons, monopoles, etc.) 
have been proposed to explain the low energy proper- 
ties of quantum chromodynamics (QCD) respecting, for 
instance, the lower part of the Dirac Operator Spectrum 
or the confinement problem 0, 0, HI 0, Q. The inter- 
est on a semi-classical understanding of confinement has 
been recently renewed by lattice studies Q. 

Other recent lattice results Q point to the semiclassi- 
cal nature of gluon correlations in the low energy regime 
of QCD. Instanton properties are traditionally measured 
on the lattice based on a geometrical localisation of in- 
stanton shapes after a cooling procedure that eliminates 
quantum fluctuations. This method incorporates a num- 
ber of known biases, such as instanton disparition and 
distortion, etc. Alternatives used by several authors are 
to fix a given number of cooling sweeps (where authors 
assume distortions are not significant) [8j or to study the 
evolution of properties with cooling (allowing to extrap- 
olate properties to the original -uncooled- situation) [9j . 

During last years, topological properties of QCD have 
been measured directly from Dirac operator spectrum, 
made possible with the use of improved Gispard- Wilson 
fermions 01 ■ O n the other hand, other method to look 
for instantonic properties that does not require of a cool- 
ing procedure arises from the study of gluon correlation 
functions in the quenched approximation, whose non per- 
turbative part can be nicely described in terms of instan- 
tons 0]. In this note, we will present some aditional 
results concerning this method, that confirm the latter 
hypothesis that gluon correlations can be described semi- 
classically. 
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Recalling some formulae from previous papers, our aim 
is to study the gluon correlations, that can be traced 
through the following euclidean scalar form factors: 
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These are the non-perturbative MOM definitions of two- 
and three-gluon Green function in Landau gauge, where 
A ^ is the gluon gauge field, T tree is the standard tree- 
level tensor for the three-gluon vertex, f a ^ a ^ stands for 
the SU(3) structure constants and the renormalization 
point is defined by k\ = k 2 = k 2 = k 2 and ki + &2 + = 
0. 

On the other hand, the gauge field in the instanton pic- 
ture and within the sum-ansatz approach, can be written 
as 

.^ = «*e^f *(^) . (3) 

where g = (6//3) 1 / 2 is the bare gauge coupling in terms 
of the lattice parameter f3, rj is known as 't Hooft symbol 
and R aa represents the color rotations embedding the 
canonical SU(2) instanton solution in the SU(3) gauge 
group, a = ].,•••, 3 (a = 8) being an SU(2) 

(SU(3)) color index. The sum is extended over all the 
instantons and anti-instantons (we should then replace 
the 't Hooft symbol rj by 77) in the classical background 
of the gauge configuration. 4>{x) is the instanton profile 
function. 

If we neglect instanton position and color correlations, 
eqs. (|2I3|) lead to 
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for m — 2, 3; where n is the instanton density. It depends 
on the functional I(kp) of the general instanton profile, 
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and < • • • > means the average over instanton sizes with 
a given normalised instanton radius distribution, p(p). 

Of course, Green functions in the instanton liquid pic- 
ture depend on the instanton profile and radius distri- 
bution and both are indeed far from being definitively 
known. As far as one searchs for some striking traces 
of a semi-classical picture success describing gluon cor- 
relations, the necessity of those two inputs disturbs our 
purposes. In 0, we solved this inconvenience not by di- 
rectly studying two- and three-gluon Green functions but 
a proper ratio of them, 
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defining the non-perturbative QCD coupling constant 
in the symmetric MOM renormalization scheme. Only 
by assuming a narrow instanton radius distribution, the 
instanton liquid model approximatively predicts a k A - 
power law for the coupling in Eq. (j^J). Lattice estimates 
of the QCD coupling constant clearly manifest such a 
fc 4 -power behaviour for the low-momentum range , al- 
though this pattern is distroyed by quantum interferences 
above k ~ 1 GeV. The latter nicely suggests that the 
larges-distances gluon correlations are dominated by the 
classical solutions of the QCD lagrangian. However, a 
rather firm confirmation of this conclusion is to be ob- 
tained after eliminating quantum fluctuations from the 
lattice gauge configurations (by a cooling procedure) and 
recovering all over the range the classical pattern. 
Moreover, the asymptotic behavior 
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is obtained with the only condition 1 of 0(0) = 1. Thus, 
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Therefore, the large-momentum limit do not depend on 
the radius distribution nor the instanton profile (the in- 
clusion of merons in the classical background would nei- 
ther modify the results in Eq. ©) and can provide a 
much cleaner indication of instanton effects. The ques- 
tion about this limit will be adressed in the following. 



1 Satisfied just by requiring a well-defined Yang-Mills action, that 
is, independent of instanton profile and even true for merons. 



III. COOLING 

After Montecarlo integration, lattice gluon field con- 
figurations arc dominated by short-distance correlations. 
No trace of instantons can be (in principle) seen except, 
as suggested by obtaining the fc 4 -power law in [||, for 
the deep-infrared domain where large-distance correla- 
tions should be dominant. A method widely used to 
detect instantons in lattice configurations is to perform 
a "cooling" procedure Q, that eliminates higher en- 
ergy modes and, as has been extensively shown, reveals 
instanton-likc structures. 

The procedure, proposed by Teper lies on an iter- 
ative change of the lattice variables by the values that lo- 
cally minimises the action. This procedure progressively 
reduces the action, suppressing quantum correlations in 
the field configurations. 

Cooling methods have been criticised as they modify 
the configuration (instanton sizes become distorted, in- 
stanton and anti-instanton anhilate to each other, ...) 
and no indication is given about how to recover the phys- 
ical information. Some gain on this problem was given by 
Ringwald et al. |l2J, that proposes that cooling acts dif- 
ferently according to the lattice spacing used. They claim 
that the effect of cooling grows as the correlations it in- 
troduces, with the square root of the number of sweeps. 
Moreover it allows to indicate the equivalent number of 
cooling sweeps for different lattice spacings in terms of 
the so called cooling radius: 
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where a is the lattice spacing and n c the number of cool- 
ing sweeps. Thus, any property of the instanton configu- 
ration can be traced through the cooling procedure as a 
function of the cooling radius and extrapolated back to 
the undistortcd configuration. 

In the next section we will use Teper procedure to sup- 
press the short-distance correlations and force the under- 
lying picture from large-distance correlations to emerge. 
We also test the evolution of the different lattice config- 
urations with cooling in terms of the cooling radius. 



IV. LATTICE RESULTS 

We have performed simulations for (3 = 5.6, 5.8 on 
a 24 4 lattice and for f3 = 6.0 on a 32 4 , obtained the 
gauge field Montecarlo configurations and computed the 
two- and three-gluon Green functions. Then, we applied 
the Teper procedure discussed in the previous section to 
"cool" the field configurations and again computed the 
Green functions for different numbers of cooling sweeps. 
We present the results in the plots of fig. ^ and, in the 
following, we will discuss them. 
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FIG. 1: Evolution o/G (2) (fc 2 ) (7e/tJ and G (s) (k 2 ) (right) with 
cooling, a clear transition happens after a few cooling sweeps. 



A. Green functions 

The first and more impressive hint we observe in the 
plots of fig. is that both gluon propagator and vertex 
join the asymptotic classical behaviour predicted by Eq. 
(|5|) after a few cooling sweeps. In particular, for the gluon 
propagator (left plots), while in perturbation theory it 
decreases at high momentum, k, with fc~ 2 (except for the 
anomalous dimension logarithm), after killing quantum 
correlations we obtain a rather clean k~ e behaviour. For 
the three-gluon vertex (right plots), correspondingly, we 
observe A: -10 instead of k~ 6 after cooling. These observa- 
tions allow to conclude that after cooling, correlations at 
high energies arc dominated by instanton-like structures. 

A second remarkable hint is how the asymptotic classi- 
cal regime is reached: the higher is the number of cooling 
sweeps, the lower are the energies where both Green func- 
tions begin to behave as given by Eq. (JSJ) . As a matter of 



fact, the typical instanton size, p, grows with the cooling 
and the asymptotic formula Eq. © stands for k ^> p^ 1 . 
Furthermore, if we compare lattice data at different (3 
for the same number of cooling sweeps, we also qualita- 
tively observe that the effect of cooling grows with the 
lattice spacing, a(/3), as expected after considering the 
cooling radius in Eq. © (the lower is /3, the lower are 
the energies where the asymptotic behaviour is reached) . 



B. Evolution with cooling 

The effect of the cooling in terms of the cooling ra- 
dius can be quantitatively analyzed because the instan- 
ton density can be extracted from a fit of the lattice 
data in the asymptotic regime to Eq. f° r au the 
lattices and number of cooling sweeps. The asymptotic 
behaviour does not depend on the details of the instanton 
liquid, and is the same for any instanton profile or even 
for merons 0. Thus, the fitted density should not be 
sensitive to the possible instanton deformations induced 
by the cooling. 

In fig. |21 we present the evolution of the density with 
the number of cooling sweeps and with the cooling radius. 
The nice matching of the results for the three lattices 
proves that the estimated density through the fit to Eq. 
(JSJ) indeed evolves with the correlation length introduced 
by the cooling. In other words, for the same cooling 
correlation length all our cooled lattice data match to 
each other and behave as classically predicted in Eq. JSJ. 
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FIG. 2: Density obtained from a fit of the lattice data for 
G^ 3 '(k 2 ) to the classical expresions Different curves for 
different values of (5 (detail) do coincide in terms of cooling 
radius. 
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C. The running coupling 

A final hint of plots in fig. to be discussed is 
why the cooling procedure, supposed to destroy only 
short-distance correlations, immediately modify the low- 
momentum range of Green functions. This is because 
of the density damping induced by the cooling. The in- 
stanton profile depends of the instanton density 0| 
and the low-momentum Green functions are affected by 
instanton profile effects diverging as k — > in the 
zero-density limit. 

A simple manner of testing this explanation, with no 
detailed discussion about the instanton profile, is ana- 
lyzing the QCD running coupling in Eq. (jHJ). As we 
above discussed, if the radius-distribution width is small 
enough, the instanton-profile effects from two- and three- 
gluon Green function very approximatively cancel to each 
other in Eq. JfjJ) and we recover a fc 4 -power law. We can 
write: 



In fig. |21 the observance of the A: 4 -power law by the 
running coupling defined in Eq. JSJ) for both low and 
large-momentum range is manifest. We present, as an 
example, the results from the simulation at (3 = 5.6 in a 
24 4 lattice, but the same is obtained from the other sim- 
ulations. For the three lattices, the shifts of the straight 
line fitting the low-momentum data respect to the one 
fitting in large momenta correspond to c ~ 1.5 — 2.0 de- 
fined in Eq. |JTU)|. in good agreement with results in |(|. 
Consequently, densities plotted in fig. [21 obtained by fit- 
ting lattice data to the asymptotic formula Eq. (JSJl for 
large momenta, correspond in practice to the same fitted 
for low momenta in |6j. 
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The low-momentum limit is taken from 0, where c ~ 1.6 
and e ~ 0.1 - 0.2. Both were estimated for 0.4 < k < 0.9 
GeV with a instanton radius distribution taken from [fH ] 
centered at ~p = 1.5 GeV and are rather independent of 
the choice for the particular profile function. 
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FIG. 3: Log-log plot of the running coupling defined in Eq. 

in terms of momenta. The two limits predicted by Eq. HUti 

are shown by straight and dashed lines respectively. 



V. CONCLUSIONS 



In conclusion, after suppression of the short-distance 
correlations by cooling, the dominance of the semi- 
classical large-distance correlations leads to a nice de- 
scription of the asymptotic behaviour of gluon Green 
functions within the instanton picture. In the low mo- 
mentum range (fc < 1 GeV), the cooling procedure only 
modifies details of the Instanton liquid like density or in- 
stanton profile. This can be seen by analyzing the QCD 
running coupling defined by Eq. 10 which shows the 
same classical pattern for both low and large momenta, as 
predicted by Eq. (|10|) . Consequently, the low- momentum 
gluon Green functions from Montecarlo gauge configu- 
rations are determined by the classical correlations, al- 
though their description within the Instanton picture re- 
quires to know the details of the Instanton liquid. 
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